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> ! Abstract 

m 

The magnetic moment of a single impurity atom in a finite free electron gas is studied in a 

^H \ 

combined x-ray magnetic circular dichroism spectroscopy and density functional theory study of 
size-selected free chromium-doped gold clusters. The observed size-dependence of the local mag- 

m 

netic moment can essentially be understood in terms of the Anderson impurity model. Electronic 
shell closure in the host metal minimizes the interaction of localized impurity states with the 
confined free electron gas and preserves the full magnetic moment of 5 u& in Cr AuJ and Cr Au, 



clusters. Even for open-shell species, large local moments are observed that scale with the energy 
gap of the gold cluster. This indicates that an energy gap in the free electron gas generally stabilizes 
the local magnetic moment of the impurity. 

PACS numbers: 75.75.-c, 75.20.Hr, 78.20.Ls, 36.40.Cg 



The interaction of localized impurity states with a free electron eras til leads to com- 
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plex phenomena such as Friedel oscillations [2| or the Kondo effect [3J. The properties of 
magnetic impurities in non-magnetic bulk metals |4| have therefore been subject of intense 
research over the last 50 years. Considerable advances in the understanding of these many- 
body effects have been made by applying photoemission, x-ray magnetic circular dichroism 



(XMCD), and scanning tunne 
or well defined quantum dots 



mg 
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spectroscopy to adatoms [5KL0J]. clusters on surfaces [11] . 



131 ] . This allows the study of basic parameters such as the 



on-site Coulomb repulsion or the amount of interaction of the impurity atom with the host 
metal. In all these cases, however, the impurity atom is in contact with a bulk free electron 
gas that has a continuous density of states. In contrast, the study of single impurities in 
size-selected clusters would allow to characterize the interaction of localized electronic states 
with a finite electron gas, i.e., with a well defined number of electrons occupying a highly dis- 
crete density of states. Because this introduces a new parameter for the control of electronic 
and magnetic properties, isolated doped coinage-metal clusters have been studied intensively 
by density functional theory (DFT) calculations [14J, Il5] and experiment: Their electronic 
and geometrical structure as well as their relative stability have been probed by photoelec- 
tron spectroscopy 16r|l8| , infrared dissociation |l9l . |20| and ultraviolet photofragmentation 
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22), as well as electron diffraction 23(. Yet, none of these experimental techniques di- 



rectly addresses the magnetic properties. Here, we use chromium-doped gold clusters as 
model systems that combine considerable local magnetic moments, carried by the 3d elec- 
trons of the impurity atom, with a finite free electron gas formed by the gold host. We study 
the magnetic moment of a single chromium impurity by local and element-specific XMCD 
spectroscopy of size-selected gas phase clusters [24|, |25|. The results of our combined exper- 
imental and theoretical study show that the size-dependent change of the local magnetic 
moment in the finite, isolated system CrAu^ is correlated with the energy gap at the Fermi 
energy of the host cluster and can essentially be understood within the Anderson impurity 
model [1]. 

Experimentally, x-ray absorption and XMCD spectra of free size-selected doped gold clusters 
were taken in ion yield mode with a combined linear ion trap and superconducting solenoid 
setup 
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27] with sufficient sensitivity to study singly doped size-selected clusters [28|, [29] . 



The cluster beam is produced in a magnetron gas aggregation source by co-sputtering of 
chromium and gold targets, mass selected in a quadrupole mass filter, and transferred into 



the radio frequency ion trap by electrostatic and radio-frequency ion guides. Inside the 
ion trap, the doped gold clusters are magnetized by an external 5T magnetic field under 
continuous helium buffer gas cooling to an ion temperature of (20 ± 5) K. The density of 
the dilute gas phase sample is ~ 5 x 10 7 ions cm -3 . Resonant photoexcitation at the L 2j 3 
absorption edges of chromium was performed at BESSY II beamline UE52-SGM. After x- 
ray absorption, the core-excited cluster ions relax via cascading Auger decay, which leads to 
photofragmentation of the parent ion. These photoion yield curves were taken for linear as 
well as for circular polarization with parallel and antiparallel alignment of photon helicity 
and magnetic field. The resulting x-ray absorption and XMCD spectra, shown in Fig. [TJ 
were normalized to the incident photon flux and scaled to unity at the L 3 edge for com- 
parison. In addition, the XMCD spectra were normalized to the number of unoccupied 3d 
states as inferred from the integrated x-ray absorption signal. These x-ray absorption and 
XMCD spectra allow us to obtain local and element specific information on the electronic 
and magnetic structure of the chromium impurity in CrAu^ clusters. This holds true even 



though the XMCD spin sum rule 30[ is not applicable to chromium because of the inter- 
mixing of 2p 3 /2 — > 3d and 2pi/ 2 — > 3d transitions. Nevertheless, quantitative information on 
the local magnetic moments can be obtained from a fingerprint analysis by comparison of 
calculated and experimental signatures in x-ray absorption and XMCD spectra. 
To analyze the electronic and magnetic properties of CrAu^, established ground state ge- 
ometries [14j] were re-optimized in a Kohn-Sham DFT framework as implemented in the 



quantum ESPRESSO 5.0 plane wave code 311], employing the Perdew-Burke-Ernzerhof ap 



proximation to the exchange-correlation functional 32J. The kinetic energy cut-off for wave 
functions and charge density were set to 680 eV and 2720 eV, respectively. Scalar relativistic 
effects were taken into account by using an ultra-soft pseudopotential of the Vanderbilt type 



33j . In addition, atomic Hartree-Fock calculations of x-ray absorption and XMCD spectra 



for an isolated Cr + ion with [Ar] 3d 5 initial and 2p 5 3d 6 final state configurations were per- 



formed in the missing program package 34 ] 



As can be seen in Fig. CD, all CrAu^ clusters under study exhibit local and total magnetic mo- 
ments, which are aligned in the presence of the external magnetic field, and therefore lead to 
nonvanishing XMCD signals independent of the detailed geometric and electronic structure. 
Furthermore, there is very good agreement of the calculated atomic Cr + x-ray absorption 



spectrum 
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361 ] with the experimental spectra of CrAu^" and, to a lesser extent, CrAu^, 
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FIG. 1. Left: Linear x-ray absorption spectra of CrAu^ clusters, normalized to the L3 maximum 
intensity. Right: XMCD spectra of CrAu^, n = 2 — 6, normalized to the number of unoccupied 
3d states. In both cases, the spectra of CrAu^~ and CrAu^~ are overlaid with theoretical spectra of 



isolated Cr^ 1 



dashed line) from a Hartree-Fock calculation. Also shown are relaxed ground state 



structures 14| (dark atom: chromium; light atoms: gold). 



indicating atomic-like chromium 3d 5 configurations in both cases. This strong localization 
of the 3d electrons at the chromium site is also reflected in screening that shifts the CrAu^~ 6 
L 3 absorption line by « 0.5 eV to lower excitation energy [37J, |38| as compared to CrAu^_ 5 . 
The XMCD spectra of CrAu,^ corroborate this interpretation: Again, the calculated XMCD 
spectrum of ionic Cr + shows excellent agreement with that of CrAu^, and good agreement 
with the XMCD spectrum of CrAu^ within the experimental signal-to-noise ratio. For 
CrAu^", a fit of the calculated XMCD asymmetry to the experimental spectrum yields an 
alignment of 0.42, which corresponds to the Brillouin value for a total spin S = 5/2 at the 
experimental conditions of B = 5T and T = (20 ± 5) K. From this agreement, a localized, 
atomic-like 3d 5 electron configuration and thus a local spin magnetic moment of 5 //b can be 
deduced for the single chromium impurity in both clusters. This is also supported by DFT 
calculations as can be seen in Fig. 0(a): In accordance with previous theoretical results [14 1 
we find local chromium spin magnetic moments of 4.85 /zb for CrAu^ and 4.75 hb for CrAu^ 
by projecting the Kohn-Sham orbitals onto atomic wave functions. 

The situation is different for the experimental spectra of CrAu^__ 5 in Fig. (TJ where no agree- 
ment with atomic Hartree-Fock calculations can be found. This suggests that the chromium 
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FIG. 2. (a) Calculated local chromium spin magnetic moments us of CrAu^ from a DFT popu- 
lation analysis (black triangles) and as predicted in the Anderson model (red bullets); (b) Cr-Au 
interaction energy E- mt ; (c) Anderson criterion (Uq + 4 J) /T; (d) energy gap AE of the gold host. 



3d electrons hybridize with gold 5d/6s states. Apparently, this hybridization leads to a 

reduction but not to a complete quenching of the local spin moment. 

This particular behavior of CrAu^ and CrAu^j" can be understood from the size-dependence 



of the chromium-gold interaction. The go 
structurally close to pure Au 2 and Au 6 [39 



d subunits of both clusters, depicted in Fig. CQ are 
, i.e., they remain nearly undistorted when adding 
the chromium impurity. In CrAu^_ 5 , in contrast, the number of Cr-Au bonds is maximized, 
and the gold host is strained and deformed in comparison to its isolated, relaxed counterpart. 
This enhanced stability of Au 2 and Au6 stems from shell closure for two and six delocalized 



6s electrons 



39[, which, in spite of a strong spd 



lybridization 40||, form a free electron 



gas confined in a two dimensional potential well 41] . Therefore, Au2,6 are known to feature 
large second differences in binding energy and wide energy gaps AE of w 2 eV at the Fermi 



energy 



Consequently, the chromium cation can be expected to interact more weakly 



in these two cases than with the open-shell gold clusters. This is indeed true, as can be 
inferred from the chromium- gold interaction energy E- int = E (Cr + ) + E (Au n ) — E (CrAu^) . 
To extract the contribution of the chromium interaction with the gold cluster from the total 
energy, E (Au n ) is calculated in the same geometric configuration of Au n as in CrAu^. As 
expected, the weakest impurity-host interactions of 2.5 eV and 3.3 eV are indeed found for 
CrAu^" and CrAu^, respectively. In CrAu^~_ 5 , E- int increases to 3.5-5 eV and indicates an 
increasing amount of covalent or metallic bonding, as can be seen in Fig. [2] (b). In contrast 
to electronic shell closure, the coordination of the impurity atom only plays a minor role 
in CrAu^ even though it does lead to a slight decrease of the spin magnetic moment with 
increasing coordination number for two-dimensional clusters with n = 2 — 5, where the co- 
ordination number is equal to the number of gold atoms and the average chromium-gold 
distances of « 2.7 A are comparable. A similar effect is observed for larger (n = 6 — 9) 
three-dimensional clusters. 
The interaction between the magnetic impurity and the free-electron host states is a crucial 

n 

ingredient for the impurity magnetic moment in the Anderson impurity model [II] . Within 
this model, the size of the magnetic moment of the impurity atom sensitively depends on 
the interplay of the on-site Coulomb repulsion, i.e., the direct Coulomb interaction of two 
electrons in the same localized orbital, and on the width 2T of the localized state. This 
width is determined by the coupling strength to the free electron states and by their density 
in the vicinity of the energy of the impurity state (l|. Since the presence of an energy gap 
AE strongly affects this density of states and thus the amount of hybridization, a relation 
between AE and the impurity magnetic moment can be anticipated. In the absence of 
interaction with the free electron gas, the impurity states E and E + Uq are separated by 
the bare Coulomb interaction Uq that preserves the local magnetic moment if Uq pushes 
the state E + Uq above the Fermi level. In the presence of interaction, virtual states are 
formed at energies E + Uq ■ n_ and E + Uq ■ n+, where n± are the occupation numbers 
of the impurity atom majority and minority states. The separation of the virtual levels is 
reduced to an effective value U e g = Uq (n+ — nJ) by hybridization of the localized impurity 
states with free electron gas states, which causes a broadening of the virtual states |l| . Here, 
(n + — n_) is the spin polarization of the localized state, which depends on the ratio of Uq/Y 
for a single-orbital impurity. A generalization of the Anderson impurity model for a five-fold 
degenerate impurity state, as for 3d elements, predicts a transition from a magnetic to a 



non-magnetic impurity state for (Uq + 4J)/T < n, where J is the intra-atomic d — d exchange 
42|. 

For a quantitative analysis of (Uq + 4J)/T, we calculated U = (U — J) (n + — n_) in a 
self-consistent scheme for the multi-orbital systems CrAu^ 43l. |44| . Here, U is derived from 
the slope of the linear response of the occupation number of the 3d impurity states to a 
rigid potential shift introduced at the impurity site 43(. Since U is obtained from ab initio 



calculations, the Coulomb interaction among the impurity 3d electrons, and all screening 
and hybridization effects are intrinsically accounted for. Hence, Uq can be determined from 
Uq = U / (n+ — nJ) + J. Since the atomic d — d exchange interaction J is only weakly 
screened, J can be assumed to be independent of the cluster size and is of the order of 



0.5 - leV 



34 



45j . The half width V of the localized impurity states is estimated as the 
weighted standard deviation of the ci-projected density of states of the impurity. The result- 
ing values of (Uq + 4J)/T obtained with constant J = (0.75 ± 0.25) eV are shown in Fig. 
|2] (c). As can be seen, the Anderson criterion for a magnetic impurity state is well satis- 
fied throughout the whole size range, which is in perfect agreement with the non-vanishing 
XMCD signal presented in Fig. CD for all CrAu^ clusters studied here. In particular, local 
maxima of (Uq + 4J)/T are found for CrAu^ and CrAu^, where the width 2T is reduced 
because of the large energy gap AE in the gold host. 
The impurity spin magnetic moment that is predicted within the Anderson impurity model 



for these parameters Uq, J, and T can be obtained by solving 

+ 0.5 (1) 



1 
n± = — arctan 

7T 



C7 » + 4J („ =F -0.5) 



r 

[l|, |42| for the values of (U + 4 J)/T given in Fig. |5] (c). As can be seen in Fig. |5] (a), the size 
dependence of the spin magnetic moment obtained from a DFT population analysis shows 
qualitative agreement with the spin magnetic moment deduced within the Anderson impu- 
rity model, and both are in accordance with the experimental data. Hence, the magnetic 
moments observed in CrAu^ clusters can essentially be explained in terms of the Anderson 
impurity model. Interestingly, however, it seems that the Anderson impurity model tends 
to underestimate the impurity magnetic moment. Beyond the simplified picture given here, 
a fully quantitative description might be obtained by an extension of the Anderson impurity 
model to explicitly include a discrete density of states that is absent in the bulk free-electron 
gas. 



Fig. [2] shows that the variation of the impurity magnetic moment is only of the order of 
10 % for CrAu^ even though the energy gap AE and the interaction Ei nt vary by a factor of 
two. This is because even in bulk gold, i.e., in the absence of an energy gap, the chromium 



impurity is magnetic 46( and carries a spin magnetic moment of 3.61 /ie [47]. The effect 
of reduced hybridization on this already large magnetic moment is limited, the more so as 
relation ([1]) does not translate the energy gap linearly to the spin polarization |l|, |42j. A 
much larger effect of the energy gap or a discrete density of states on the local magnetic 
moment should be observed in systems that are closer to or even below the threshold of 
quenching in the bulk case. In this case, the energy gap might even serve to restore the local 
magnetic moment of the impurity atom. 

In summary, the experimentally observed size dependence of the XMCD spectra of size- 
selected chromium-doped gold clusters are in line with the Anderson impurity model. The 
size-dependent variation of the spin magnetic moment can be linked to the amount of hy- 
bridization of the impurity with the host density of states and is governed by the energy gap 
of the host gold cluster. Electronic shell closure in the gold host leads to wide energy gaps 
AE in the free-electron states, which reduces the interaction with the impurity and causes 
the maximum spin magnetic moments of 5 /lib for CrAu^ 6 . This effect is a result of quantum 
confinement and is unique to finite systems. 
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